Gravity theories, Transverse Doppler and Gravitational Redshifts in Galaxy Clusters 
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There is growing interest in testing alternative gravity theories using the subtle Gravitational 
Redshifts in clusters of galaxies. However, current models all neglect a Transverse Doppler redshift 
of similar magnitude, and some models are not self-consistent. An equilibrium model would fix the 
Gravitational and Transverse Doppler velocity shifts to be about 6<r 2 /c and Za 2 /2c in order to fit the 
observed velocity dispersion a self-consistently. This result is from the Virial Theorem for a spherical 
isotropic cluster, and is insensitive to the theory of gravity. In any case, a gravitational redshift 
signal cannot directly distinguish between the Einsteinian and f(R) gravity theories, because the 
mass of the cluster dark halo must be treated as an unknown fitting parameter, whose value must 
vary according to the theory adopted, otherwise the system would be in equilibrium in one gravity 
theory and out of equilibrium in another. 
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The theory of gravity has been subjected to various tests with the ever-improving quality of data from cosmology, 
galaxy clusters, galaxies and the solar system[l[. As shown by recent numerical N-body simulations on f(R) type or 
scale-coupled gravities Q, dynamical data on non- linear cluster scales help to break theoretical degeneracies on linear 
cosmological scales, and overcome statistical uncertainties in observations. Past techniques often propose comparing 
lensing data and kinematic data with simulations which often involve significant amount of modeling of the mass 
distribution before setting indirect constraints on the gravitational potential $(r) of the cluster. It would clearly 
be better to measure the gravitational potential in a galaxy cluster directly and compare this potential with the 
prediction from the Poisson equation for the mass distribution in a given gravity theory. 

Indeed the gravitational potential is an observable. For example, spectral lines from the surface of the Sun are 
observed to be gravitationally redshifted by a small amount GMq/ Rqc ~ 0.6 km s -1 , and the shift increases for more 
compact stars. On cosmic scales, the deepest potential well $(r) is in a cluster of galaxies with usually a bright 
central galaxy (BCG) near the bottom of the potential. A nearly spherical distribution of many hundreds of galaxies 
orbit around the BCG, with a Gaussian dispersion of random velocity of a(r) ~ lOOOkms -1 in each direction. The 
observed line-of-sight Doppler shifts of galaxies relative to the BCG satisfy a Gaussian distribution with a small but 
non-zero mean velocity. This is partly due to the Gravitational Redshift (GR), a feature in any metric theory of 
gravity, caused by the spatial variation of the gravitational potential: 



[<J>bcg ~ Hr)] fc (1) 



This signal of ~ lOkms -1 becomes detectable above the a/\/N uncertainty of the mean velocity once the sample 
size TV > 10 4 . To obtain such a large sample for the first time, Wojtak et al. 0| use TV ~ 78000 galaxies from about 
eight thousand clusters from the Sloan Digital Sky Survey (SDSS), divide the galaxies into four bins according to 
their projected distances R from their respective BCGs, "stack" their light-of-sight velocities relative to their BCGs, 
carefully remove interlopers, and finally compute the mean velocity in each bin. In this Letter we investigate the pros 
and cons of the gravitational redshift approach, and for the first time introduce a new effect in galaxy clusters. 

In fact, the Gravitational Redshift is supplemented by an additional redshift of comparable amplitude. For any 
metric theory of gravity [H the space time near a galaxy cluster is described by the metric dr 2 = (1 + 2Q/c 2 )dt 2 — 
(1 + z)~ 2 (l + 2^ / c 2 )dx 2 . Light emitted from a cluster at redshift z is time-dilated so that the ratio of the observed 
wavelength to the emitted wavelength changes as: 



±*S- = il + z) 



v 2 /2 
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so we have an additional shift 



[<|v| 2 )-|vbcg| 2 ]/2 C , (3) 



owing to the Transverse Doppler (TD) effect from random motions of galaxies in special relativity (SR) . Here we have 
neglected the much smaller shifts caused by the potential and motions within individual galaxies or stars. Wojtak et 
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al. [J reported a blueshifting of the mean apparent line-of-sight velocity of the galaxies in the SDSS clusters, again 
relative to the BCG, which was then interpreted as GR. But this interpretation is incomplete. The TD effect always 
co-exists in proportion to GR because of the Virial Theorem: 

<-$/2) /2c ={GM/r) /2c = (|v| 2 ) /2c (4) 

where () denotes the averaging over all gravitational masses in the whole virialized volume of a cluster, and the factor 
of 1/2 in front of $ prevents double counting of the pairwise mutual potential. Thus the random kinetic energy 
per unit mass v 2 /2 is globally 25% of the average potential — <f>. The ratio of 1/4 holds even after averaging over a 
distribution of clusters of different mass and for clusters of any density profile and anisotropy parameter, so the Virial 
Theorem is a robust link between TD and GR effects, and their superposition is observed as the mean velocity shift. 

Recent tests of general relativity are often posed in the context of the A Cold Dark Matter (CDM) model. It 
seems straightforward to test many gravity theories with their gravitational redshift prediction, but the invoking of 
dark haloes complicates the comparison of some gravity theories. E.g., it is well-known that in f(R) gravity Q the 
Newton constant G can be boosted by a factor which varies between and 4/3, depending on the environmental 
matter density. For a fixed cluster mass Wojtak et al. claimed that the 4/3 ~ 1.33 boost of the Gravitational Redshift 
signal in the Hu & Sawicki f(R) model Q with |/ro| = 10~ 4 robustly distinguishes Modified Gravity from Einsteinian 
Gravity. Such a claim, however, has a flaw: the total mass of the dark halo is an unknown free parameter, which must 
be determined by fitting the observed velocity dispersion as a function of distance from the cluster centre. Since G 
appears only in the combination GM, one can cancel essentially the enhancement of G in f(R) gravity by reducing 
the halo mass parameter M, thus obtaining the indistinguisble fit to the velocity dispersion curve and to the mean 
velocity shift signal. 

A more specific example is to use the hydrostatic equation 



d<£ d(na 2 ) 
dr ndr 



(5) 



For an isotropic number density of tracers (here galaxies) n oc r 1 at large radius r for a Keplerian potential 
& = -GM/r. This yields the solution 

2c a [T) ~ 2c r J r r 2 r~< ~ 2( 7 + 1) rc ' [ ' 



which implies that the random kinetic energy v 2 /2 = 3er 2 is locally ~ 3/8 — 3/10 of a Keplerian potential | — $(r)| 
of a galaxy count profile with 7 ~ 3 — 4 at large radii. The ratio 3/8 or 3/10 holds even after stacking of clusters 
of different masses and line-of-sight projection. This argument is true in standard gravity as well as in f(R) gravity. 
Clearly gravity models with the same GM predict the same dispersion curve, and velocity shifts. 

To compute the TD and GR effects generally at any projected radius, we start with the isotropic Jeans equation 

d(na 2 ) d<S> 

for the observable tracers (galaxies) with a number density n(r) in equilibrium in the potential $(r). We integrate this 
equation over the line-of-sight depth Z through a cluster after multiplying by ZdZ on both side. Using the integration 
by parts, Zd{na 2 ) = —dZ(na) + d(Zno 2 ), where the total derivative term vanishes at boundary, we find 

dZ {no 2 ) = J dZ [nZ^j . (8) 

Multiplying by 2nRdR on both sides of the above equation predicts that the specific 3-dimensional kinetic energy 
averaged in an aperture of projected radius R to R + dR is (|v| 2 ) /2 = (QZ(d<fr/dZ)) = (GM(r)QZ 2 r~ 3 ) inside the 
virial radius. Here Q = 3/2 from summing the three velocity components in quadrature. This equation allows us to 
predict the SR effect at all radii for any matter density in any metric-based gravity theory, since the Jeans equation 
applies to any force which is a gradient of a potential. E.g., if the density ~ r~ 7 ~ r -3 and gravity ~ r~ 2 then 

= SttGM 

2 8 4R { ' 



Here the factor irGM/AR comes from the density-weighted line-of-sight integration of —$(r). 
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To compute the GR and TD effects for the SDSS clusters, we follow the same methodology as Wojtak et al. We 
account for different cluster masses using a Salpeter-like mass function dN/dM vn - ~ M vir ~ M vil taken between 
the mass range Mi and M u . 



(|v| 2 )/2 = A(R)- 1 / dM vil / dZnWfQ^r-^d^/drJjl^^gs+sr, (10) 
($) = A(R)- 1 dM vix ^n(r)$(r)| r=V23W , (11) 



where = J* M " dM v j r dZn(\/ Z 2 + R 2 ) is essentially a stacked line-of-sight integrated density at the projected 

radius i?. and the spherical potential and tracer (galaxy) number count density are given by 

CM ■ 

$(r) = ln(1 + r C/r vil ), (12) 



where Ci 



ln(l + C) - 



n(r ) BM ^ /3 ^ fl3 ) 

() ~ AirCi r(r + r vir /C) 2 ' 

We fix the halo concentration parameter C = 5 and the virial radius r v - n - = 



1 /3 

1.2 (M v i r /10 14 M Q ) , as in Wojtak et al. Such a spherical NFW potential is an approximation to the true po- 
tential in Einsteinian gravity; in f(R) gravity the potential from N-body simulations tends to be more concentrated 
Q and in TeVeS gravity the potential tends to have a pure ln(r) profile at large radii. For the lack of a more precise 
analytical model of the tracer (galaxy) count n(r), we have used the spherical NFW profile, and assumed the number 
of observable galaxies inside the virial radius to be N v - lT ~ M vir /cst. In reality, the number count profiles of observ- 
able bright galaxies tend to be more concentrated than the halo mass density. The selection criteria of galaxies of 
measurable redshift should be modeled as well, e.g., the preference for brighter and nearer galaxies and the velocity 
cut in the line of sight makes n(R,Z) cylindrical instead of spherical, and J dM v i v J^ 0Mpc 2dZ f® Mpc 2irRdRn(R, Z) 
is equal to the total size of the sample, about 78000 galaxies for Wojtak et al., which fixes the normalization con- 
stant B. Finally the BCGs are not always in the centre. We shall neglect these complexities, but empirically let 
Q ~ (3r v j r + r)/ (2r V j r + r) < 3/2 to mimic the effects of mild anisotropy and non-equilibrium at large radii (e.g. Fig. 
1-3 of i). 

The results for Einsteinian gravity gravity is shown in Fig. 1, where we assume a halo mass range of (M/,M U ) = 
(0.11 x 10 15 , 2 x 10 15 )M Q . Note that these are fitting parameters deduced from hydrostatic balancing of the pressure 
from the observable velocity dispersion a and the halo gravity, as one cannot directly observe the halo and measure 
its mass. One can see our choice of parameters can fit the observed (3(T 2 bs — 3cr BCG )/2c curve and, in doing so, we can 
predict the — |$(i?) — ^bcg)!/ *-*R curve - Note the inevitable reversal from the observed average ~ 6.5 ± 4kms _1 
blueshifting to redshifting when within 0.2 Mpc of the BCG (cf. grey error bars and lines in Fig. 1) due to TD: 
the line-of-sight dispersion of non-BCGs a bs(R) — 600 km s -1 ~ 3<7bcg converts directly to a {3a 2 bs — 3ct bcg )/2c ~ 
1.6 km s" 1 TD differential redshift near an isotropic cluster centre. The TD signal (red crosses) is clearly both 
non-negligible and model-insensitive, is thus a robust constraint applicable to any metric-based gravity theory. 

As stated earlier, one should not compare an f(R) gravity model with an Einsteinian gravity model for with 
the same halo mass distribution, namely (Mi,M u ) = (0.11 x 10 15 ,2 x 10 15 )M Q , since it would overpredict the 
velocity dispersion curve a 2 (R) everywhere by the same factor of 4/3, which can be ruled out even without measuring 
gravitational redshift. This is essentially because a unique relation between the SR and GR effect due to the isotropic 
Jeans equation. The more relevant mass distribution, where all (halo virial) mass is lowered by the same factor 4/3, 
would predict a velocity dispersion curve identical as the Einsteinian curve. Instead, to show some difference, here 
we adopt (M h M u ) = (0.09 x 10 15 ,1 x 10 15 )M o , and the result is shown as solid lines in Fig.l. This f(R) model 
produces gravitational redshift and TD shifts by an amount essentially identical to Einsteinian gravity. So GR, TD 
and the velocity dispersion profile contain essentially three redundant copies of information about a metric theory, up 
to some uncertainty of anisotropy. 

Likewise, the claimed (0 — 10)kms _1 extra shift in TeVeS reduces to only (0 — 3)kms _1 when Q adopting mass 

T 2 



models consistent with cr 2 hs . Unfortunately the Transverse Doppler effect is left out explicitly in all these papers; e.g., 

ha.l-indenendent SR red-ward shift of / 2-ZJL$>(rW = ■ 

■yc 



TeVeS predicts a roughly radial-independent SR red-ward shift of {^Z '^$(r)} = ~ 1 kms 1 for 7 ~ 3. Further 
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FIG. 1: Consistency of the data with the competing GR effect — [$(r) — $bcg]/c plus TD effect (v 2 — v| CG )/2c; the 3D 
kinetic energy v 2 /2 at any R equals QZdQ/dZ averaged over the line of sight depth z, where Q = 3/2 within the virial radius 
r v i r of an isotropic cluster; slight reduction of the TD signal (yellow line) if using a smaller Q ~ (3r v h- + r)/(2r v i r + r) to mimic 
any mild anisotropy and non-equilibrium at large radii ((H). All models use the NFW halos weighted by M~ 7 ^ 3 in the virial 
mass range M = (0.09 — 1) x 10 15 M Q in |/ho| = 10 -4 gravity (solid), where the effective gravitational constant (Fig. 3 of Q) 
in these low mass halos is boosted by 33%, or M = (0.11 — 2) x 10 Mq in Einstein gravity (dashed). 



investigation including all relativistic effects and varying halo parameters perhaps within ranges from N-body f(R) 
simulations [H, 0] is necessary before making a robust evaluation of gravity theories with this new probe. 

Finally there is also the future possibility to use the spectra of cluster gas to measure the GR and TD effects 
for every X-ray cluster. The o j \f~N uncertainty in the mean velocity will be negligible for the countless ionized gas 
particles, and there is less need for stacking different clusters. The signal will differ from the technique of Wojtak et 
al. because the X-ray gas particles have a profile different from that of galaxy number density, and less anisotropy 
than the velocity distribution of galaxies. By comparing the Transverse Doppler signals of different tracers one might 
even be able to infer the anisotropy of the motion of galaxies inside clusters. 
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